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The work is devoted to the investigations of possible observational manifestations of protoobjects related to 
"dark ages" epoch (10 < z < 1000), before formation of self-luminous galaxies and stars. These objects can distort 
the cosmic microwave background. Formation of these objects is described in the "pancake theory" and in the 
model of "hierarchic clustering" . According to these theories we may consider these protoobjects as flat layers. 

We consider both Thomson (with Rayleigh phase matrix) and resonance (for complete frequency redistribution) 
scattering of cosmic microwave background radiation by moving flat layer. The appearing anisotropy and polar- 
ization of cosmic microwave radiation are calculated for a wide range of a layer optical thicknesses (from optically 
thin layer to optically thick one). Analytical solutions are also obtained for the case of an optically thin layer and 
are compared with the numerical ones. 

1 Introduction 

We consider a special period in the Universe evolution extending from recombination epoch {z = 1000) till the time 
of self-emitting stars and galaxies formation [z = 10). This period is called "dark ages" and the question about 
possible observational manifestations of dark objects from that period is actual. One of the most effective ways 
to solve this problem is to try to observe the cosmic microwave background (CMB) distortions caused by these 
objects, as it was firstly reported in PQ. The character of such distortions is defined by the physical properties of 
"dark ages" objects, which are naturally the protoclusters of galaxies 

Now two basic theories of galaxy clusters formation are being discussed. The first one, so called "pancake 
theory", was suggested by Ya.B. Zeldovich (see 01 for example) and further developed in some papers (see 0], |S]). 
According to this theory, the substance is collapsing into flat formations — "pancakes" (protoclusters). Initial 
shock waves forming in this process can be modelled as flat layers. According to the second theory, "the model 
of hierarchic clustering" , small structures are merging into larger ones ( "the boxes" ) . These structures contain 
a large number of small dense objects merging to form present galaxies and diffuse gas substance. The difference 
between these two galaxy clusters formation models consists only in the value of mass ratio of dense objects and 
diffuse substance. In the first theory gas "pancakes" contain the most part of the total mass, and in the second 
one the most but not the whole mass is contained in dense objects. Therefore, the gas mass is not equal to zero 
and its spatial distribution as flat "pancakes" takes place in both theories and we may consider them as flat layers. 

Possible distortions of CMB appear due to the fact that these flat layers can move relatively CMB with peculiar 
velocities (see 7^ for example). Due to the Doppler effect caused by this motion the CMB becomes anisotropic 
in the reference frame of the layer. This anisotropy, which is axially symmetric with respect to the direction of 
the motion, leads to polarization and intensity changes of scattered radiation. These distortions belong to the 
"secondary CMB distortions" . Their differences from ones forming at the hydrogen recombination time (z — 1000) 
consist in other angular scales and possibly in a more complicated spectrum. In particular, "secondary fluctuations" 
can appear not only due to the scattering on electrons, but on primary molecules as well |H], [HIj [TU]. 



Peculiar velocities of substance large scale fluctuations for a wide range of cosmological models can be estimated 
by the formula [TT ] . [H ) : 

Vp = 600 {1 + z)-^/'^ km/a, (1) 

where the numerical coefficient is obtained as a result of a large number of velocity estimations of nearby clusters of 
galaxies. The velocity dependence on z is defined by the general law of the massive body motion in the expanding 
Universe. The calculation of the optical depth is different for scattering on free electrons (Thomson scattering) and 
for resonance scattering on atoms, ions or molecules. Thomson optical depth does not depend on the wavelength 
and is defined by the size of density fluctuation and electron density. The last quantity sufficiently depends on the 
redshift: after the hydrogen recombination at z = 1300 the fractional electrons density falls from the unity to the 
value about 0.01% at z = 100. But after that time the number of free electrons rises due to the secondary ionization 
the details of which are being investigated by several authors and may sufficiently differ in various models (see 
for example). Observational restrictions based on the power spectrum of primary CMB fluctuations give for the 
average optical thickness of such objects tq < 0.1. However, it can achieve unity for particular objects. 

Existence of free electrons and protons leads automatically to the appearance of bremsstrahlung not related 
with the peculiar velocity. The role of such emission in the secondary ionization epoch was investigated in jl4) 
where it was pointed out that there is a quite wide range of fluctuation parameters for which the role of this process 
is very small. 

As concerned a resonance scattering in molecular lines there is a strong dependency of observed effect on the 
frequency. The estimation of the optical thickness of the layer (averaged over the line) tq in this case is much more 
complicated problem since the density of different molecules strongly depends on physical processes in a given 
epoch jnj. In principle, the estimation of second type shocks efficiency (excitation and disactivation by collisions 
with electrons and atoms) is necessary here. However, we may neglect such processes for the considered epoch 
because of a very low absolute concentration of colliding particles (numerical estimates are made in 0). Thus, in a 
first approximation we may consider only the Doppler mechanism of CMB distortions formation. CMB interaction 
with molecules having high dipole moments will be most effective. From this point of view the main attention was 
paid to the molecules HD~^, HeH'^, LiH and some other ones ^Hj- A large number of HeH^ molecules can 
be formed in the shock wave at the stage of nonlinear collapse of primary "pancakes" . This way to observe the 
large scale distribution of substance at z ~ 10 appears to be possible. The role of other molecules is many times 
less (see ^U] for example). The lines of rotational and vibrational transitions of such molecules at z ~ 200 can be 
situated in a millimeter range. The discovery of two lines of rotational-vibrational array should allow to define the 
type of molecule and z. 

The effects of CMB scattering on free electrons in the hot intergalactic gas of rich clusters of galaxies were 
considered in the papers by Zeldovich and Sunyaev ^^l) QH|i where the analytical estimates of CMB polar- 
ization caused both by the gas motion and by the intrinsic anisotropy of CMB were obtained. Those papers were 
principally based on a single scattering assumption which, however, leads to polarization which is proportional only 
to [v/cY . Analytical estimate of 0{v/c) effects with account of double scattering for moving spherical objects was 
made in |l()j . 

The primary goal of our work is the numerical calculation of CMB distortion effects depending on the layer 
parameters with no restriction concerning to the number of scatterings taken into account. We consider both 
Thomson (with Rayleigh phase matrix) and resonance (for complete frequency redistribution) scattering of cosmic 
mircowave background by moving flat layer. The appearing anisotropy and polarization of cosmic microwave 
radiation are calculated for a wide range of a layer opticall thicknesses (from optically thin layer to optically thick 
one). Analytical solutions are also obtained for the case of an optically thin layer and are compared with the 
numerical solutions. Since the electron number density after the recombination epoch becomes too small, this 
mechanism will be more effective at the beginning of nonlinear stage of protoobjects evolution at the "secondary 
ionization" epoch (10 < z < 50) (see QT] for example). The molecular scattering can play sufficient role at z < 150. 

The predicted effects can be observed in a wide range of angular sizes which depend on cosmological model 
parameters and may vary from dozens of angular seconds to dozens of angular minutes. 



2 Rayleigh scattering 

Let us consider plane-parallel horizontally homogeneous layer moving as a whole with a velocity v (with regard to 
CMB) at the direction of an outer normal n to the layer. According to eq. peculiar velocities of such objects 
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are about 150 km/s. Thus, we may assume vjc ^ 1 and disregard the terms of the order of {v/cY in further 
calculations. The layer is illuminated by CMB with the intensity (in its own reference frame) described by Planck 
function with the temperature T = 2.7K. Having transformed to the layer reference frame the CMB intensity is 
given by 

_ 2hv^ 1 

^ ^) ~ gi{l+tif3)hu/kT _ I ' 

where v is the radiation frequency, fi is the cosine of the angle between the velocity and radiation propagation 
directions in the layer reference frame, /3 = v/c, 7 = 1/ ^yl — f3^ . 

Since we have azimuthal-symmetric picture, the field of radiation is completely described by two-component 
Stokes vector i — {I^Q)"^, where "T" means matrix transposition, / = I{r,i.i,i'), Q ~ Q{r,fi,i'), r is the distance 
from one of the layer borders. 

In the case of Rayleigh scattering we will use the optical depth r as the geometric variable instead of r: dr = 
—UeCrdr and denote the layer optical thickness as tq. Here Ue is the electron density, a is the Thomson cross-section. 
Since the picture is independent of frequency, the Stokes vector can be written as i(T, /i) = (/(t, 11), Q{t, /i))'^. This 
vector is the solution of radiative transfer equation (^1], pUj l 

^dii^ = i(r, m) - (1/2) Pifi, M')i(T, fi')dfi' ~ s* (r, ^i) ^ i(r, fi) ~ s(r, ^) (3) 



dr 

with the boundary conditions 

i(0,/i) = ii(/i), ^i<0] i(To,Ai) = i2(A*), M > 0, (4) 

where ii(/i) and i2{fJ.) are the Stokes vectors of radiation illuminating the layer from outside at the boundaries 
T = and T — Tq, respectively. Here s* characterizes primary sources distribution in the layer and the phase 
matrix P{fi, /i') is, in general, the superposition of the Rayleigh phase matrix and the phase matrix of isotropic 
scattering Pi: P = {1 — W)Pi + WP^. Here W is the depolarization parameter which is usually in the range [0,1]. 
For the phase matrix the following factorization is obtained (see [201): P{ii,ii') = A(/i)A'^(/i'), where the matrix 



Here and below b = ■^/W/S. In the case of Rayleigh scattering considered here the depolarization parameter W is 
equal to 1. 

Taking into account a single scattering of an outer radiation and using the phase matrix factorization, we reduce 
the problem to the solution of eq. © with the free term 



s*(r,M) = (l/2)i(M) / i^(M') e-(-«--)/^'i2(Ai')+e-"/'''ii(-M') dp.' (6) 
and zero boundary conditions 

i(0,/i) = 0, /i<0; i(To,/i) = 0, Ai>0. (7) 

For ii(/i) and i2(/i) we use the expansion of background radiation in the layer reference frame ^ with the accuracy 
of the order of 0{v/c): 

ii(/i) = hip) - B{iy, T) [1 - {v/c)a,y] ei, (8) 

where = xe^ / [e^ — 1), x = hv/kT, ei = (1, 0)"^. 

According to eq. (jHJ we devide the both sides of the main equation © by the depth independent factor B{v, T) 
so that below the Stokes vector i(T, /i) and the vector source function s(t, /i) are dimensionless (measured in the 
units of B(i',r)). 

The factorization of the phase matrix P(/i, /i') leads evidently to factorization of the source term in the 
righthandside of eq. s(t, ji) — A{^)s{t). Then we use the procedure described in to reduce the problem to 
the solution of the linear integral equation for a vector source function s(t) depending only on the optical depth. 



3 



1.0 



0.8 - 



0.6 - 



3 



0.4 - 



0.2 



0.0 



1 

- 


1 1 1 1 1 
a 


1 1 1 


1 




























0.25 












0,08 - 



0.06 - 



0.04 - 



0.02 



1 


1 1 1 1 1 1 1 1 1 1 1 1 1 

b 




0.1^^_^_ — " — " — — " 
0.075__ ; 


- / x' 


0.05 




0.025 




0.01 






0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



0.00 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 1: Profiles of intensity change for Rayleigh scattering for different values of tq marked near the curves. 



Finally we obtain the following expression for the Stokes vector of radiation emerging through the boundary r = 
(for the observer reference frame): 



(9) 



i(0, Ai) = ei + {vlc)a^ii (l - e"^"/^) ei + (w/c)aU(/^) j er^ ' ^'^(T)dT j ^i, ^ > 0, 
where s(t) = (si(t), SQ(r))'^ satisfies to the vector integral equation 



with the core matrix (see [201 ) 

k{T) = 



s(t) = (1/2) / K{\t ~ T'|)s(T')dT' + s*(r) 
Jo 



b[Ei{T) - 3^3 (r)] 2b^[5Ei{T) - UEsir) + QE^ir)] 



(10) 



(11) 



(12) 



and the primary source term 

s*(r) = (l/2){i?3(r) - E3{t„ - T),b[E3{T) - E^ir^ - t) - i{E^(T) - E^{t^ - r))]}^. 

Here b ~ -^/l/S, Enir) — / e^'^/^/i"^^(i/i is the n-th integral exponent. 

^0 

For the intensity change AI / Iq = (I — B) / B and for the polarization degree of radiation emerging through the 
boundary r = we obtain from eq. (j^ 



A///0 = {v/c)a,T{ii,T^) + 0{{v/cY), 
Q/I = -{v/c)a^Pi^i, To) + 0((i;/c)2), 

where the profiles of intensity change and polarization are 

rifi, To) - f 1 - e-^"/^) + r e-^l^ [si(t) + (l/\/8)(l - 3m')sq(t)1 dr/^i, 



(13) 



P(m,to) =-(3/V8)(1-m2) / e--/^SQ(T)dr/M. 

Jo 



(14) 
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Figure 2: Profiles of polarization for Rayleigh scattering for different values of tq marked near the curves (a) or 
listed in the same order as the curves follow, from up to down, (b). 



For the maximum polarization degree which is reached at = and for the corresponding intensity change we 
have from eq. (|14|l 

P(0,to) = -(3/V8)sq(0), 
r(0,ro)=si(0) + (l/V8)sQ(0). ^'^^ 

For the radiation emerging through the boundary t — tq we obviously should change the sign of v in the 
formulae above i.e. to change the sign of polarization and intensity variation. 

Integral equation (|10|l was solved numerically by discretization at some r grid and reducing to the system of 
linear algebraic equations for the values of the source function at the knots of the grid. These values are used then 
to calculate the profiles of the intensity change and polarization through eq. (|14() . The results are shown in Figs. 
1 and 2. 

As we can see in Fig. 2b there is a strong dependency of polarization degree on fi for small tq. For tq <C 1 
we obtain the following expansions for si(t) and sq(t) with the accuracies up to the first and second order on the 
optical depth respectively: 

( ) - f ^ f -T + To/2 \ 

V y V (To/8V2)[(To-2T)(C + l/2)-Tlnr+(To-r)ln(ro-T)] J' ^'^^ 

where C = 0.577216 is the Euler constant. Substitution of these expansions into eqs. (|14|l and IjlSfl gives us the 
analytical expressions of the intensity change and polarization profiles for tq <C 1: 

, . ( To/2, fl < To, 

[To, ^ » To, 

P(n.) r P(0,ro)^-(3/32)r2(lnTo + C+l/2), /i « ^"0, 

\ -(l/64)r4(l-^2)^-2(i„^^_^.C. + 2/3), /i » ^-q. ' ^ 

Eq. (|18|l explains the strong dependency of polarization degree on /j. It is confirmed by numerical calculations. 
For example, at tq = 0.01 its error is less than 4% for /i > 0.09, at tq = 0.025 is less than 4% for fi > 0.14 and at 
To = 0.1 is less than 10% for /x > 0.13. We can obtain even more precise approximation taking into account the 
terms of larger order with respect to to. Thus, for P(0,to) we have 

P(0,To) «-(3/32)T2[(l + 1.25To)lnTo + C+l/2], (19) 



5 



the error of which is less than 1% for tq < 0.1. 

An approximate analytical estimate of the maximum polarization degree for a homogeneous sphere with a 
small optical radius tq was obtained by Sunyaev and Zel'dovich ^^1^ (Q/-f)max — ± (wt /c)ajy Tq /lO, where vt is 
the transversal velocity. However, the accuracy of this estimate is unknown because it was not compared with 
numerical calculations. 

Obtained dependencies of CMB intensity and polarization on layer parameters allow us to model observational 
manifestations of these objects. In the protoobject directions we would registrate increase or decrease of CMB 
intensity depending on the sign of velocity v. Having measured the intensity change, we may estimate the {v/c)to 
value through the eqs. H13|l and (|17|l . Numerical analysis of polarization in couple with intensity variation data 
would permit to obtain the values of velocity and optical thickness of the layer separately. 



3 Resonance scattering 

Let us consider now the case of CMBR scattering in a spectral line, say, in a line of HeH~^ with the laboratory 
wavelength 149/z. At z = 150 this corresponds to the wavelength 2.25 cm. 

We use the model of two-level atoms with the transition frequency 1^12 ■ Let x be the dimensionless frequency: 
X — {v ~ 1^12)/ Ai^Oi where Avo is the Doppler width of the line. Let (pix) be the absorption coefficient profile 
normalized as follows: 4>{x)dx = 1. 

Since we have axial symmetry again, the Stokes vector i will have only two components (/ and Q). But now it 
will depend not only on coordinate and direction (r and ^) but on frequency x. The radiative transfer equation is 
written as (see |OT for example): 

^^}!^Z^^2^ ^ ij)[x)i{T,ijL,x) [ d^' [ R{n,x]fi',x')i{T,n',x')dx'-s*{T,fx,x) = 



= (j){x)i{T,^i,x) - s{T,ii,x). (20) 

The matrix R(fi,x; fj,' ,x') describes the redistribution over frequencies, angles and polarization conditions at a 
single scattering, A is the single scattering albedo, s*{t,^,x) is the primary source function vector. If there is no 
primary sources embeded in the layer it is defined by illumination from outside: 

A f°° r , , 1 

s*{T,fi,x)^- dx' i?(^,a;;Ai',a;') e-("°-")"^(")/^i2(M',:E') + e-"*(")/^ii(-Ai',2;') V, (21) 
2 J~oo Jo '- 

where ili^^,x) = i2i^^,x) ^ B{vi2,T)[l~ (t;/c)a^i2^]ei, a^.^^ = ^12/(1 - e^^^^), 0:12 = hvi2/kT, ei = (1,0)"^. As in 
the case of Rayleigh scattering we use below dimensionless Stokes vector and vector source function (measured in 
the units of B{vi2,T)). 

We use the assumption of complete frequency redistribution, according to which (see |21j ) R(^,x; ^' ,x') — 
<j){x)(t){x')P{^, /i') = 4){x)4){x')A{^)A^ [jj!), where A{^) is defined by eq. (0). The value of depolarization parameter 
W in that equation is defined by the quantum numbers of the total angular momentum of the upper and the lower 
levels of transition (see [T^ for example). Polarization becomes smaller as W decreases. W = \ corresponds to 
dipole scattering. In this case the phase matrix is the same as for the Rayleigh scattering. 

The factorization of the redistribution matrix R{ii, x; fi' , x') leads to factorization of the source term in the 
righthandside of eq. (|20|l : s{t,ij,,x) — (P{x)A{ij,)s{t). Then we use the procedure described in |2] to reduce the 
problem to the solution of the linear integral equation of the type H1U|I for a vector source function s(t) depending 
only on the optical depth. But now the elements of the core matrix K{t) in eq. H10|l are (see (21)'): 

/oo 
cl,''{x)Ei{^{x)T)dx, 

a2 



Ki2{r) - K2i{t) - bX / [Ei{^{x)r) - 3Es{^{x)t)] dx, (22) 

J — CXj 

K22{t) - 262a / ^"^{x) [5Ei{4>{x)t) ~ 12E3{^{x)t) + 9E5{4>{x)t)] dx. 
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Figure 3: Angular profiles of polarization in a center of a spectral line for different values of tq marked near the 
curves (a) or listed in the same order as the curves follow, from up to down, (b). 



Finally we obtain the following expression for the Stokes vector of emerging radiation: 



iiO,^i,x) = ei + {v/c)a^,^n (l - e ei - (1 - A)[l + {v/c)a^,^^l]A{p) j e~^/«So(r)dr/C+ 



+ {v/c)a,,,A{^i) / e-"/?s(r)dr/C, M > 0, 



(23) 



where C — fi/<j>{x), so(t) is the solution of type (|10|l equation with the free term s* = ei, and s(t) = (si(t), sq{t))'^ 
is the solution of the same equation but with the free term s* = (sj , Sq)"^: 



(x) [E3{^{x)t) - E3{cj){x){To - t))] dx, 



(24) 



{x){[E3{^{x)t) - E^{c^{x){to - t))] - 3 [E^{<I>{x)t) - E,-,{<j>{x){To - r))]} dx. 



(25) 



Since for the most part of molecules the value of 1 — A is about 10"^ and the value of v/c is about 10^'', we 
may neglect the term proportional to 1 — A in eq. H23I) . However, the account of this term does not bring any 
difficulties. 

As in the case of Rayleigh scatering we introduce the intensity change and polarization profiles of the emerging 
radiation at r = in analogy to eq. (|13|l where now according to eq. (|23|) 



r{x,n,To) = n 



\ _ g-0(2:)-ro/M 


/■To 

+ 0(x) / e- 


4>(x)T/ti 




"'0 





Si 



(r) + vnl78(l-3/i2)sQ(T) dT|^l, 



P{x, To) = -3^/W/8 (1 - fi^)(t,{x) / e-*(^)^/''sQ(T)dT//x 



(26) 



For the maximum polarization which is reached at /i = and for the corresponding intensity change we have 
from above equations 

F(x,0,ro) = -3^/f sq(0), 

(Zl) 



r(x,0,ro) = si(0) + JfsQ(O). 



Note that these quantities do not depend on x. 
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The integral equation with the core matrix and free term, given by eqs. (|22() and (|24I25() respectively, was 



-1/2, 



and 



solved numerically as in the case of Rayleigh scattering. We have used the Doppler profile (j){x) - 
brought W — 1, X = 1. These values of parameters W and A evidently maximize the degree of polarization. Fig. 
3 shows the dependency of polarization P on the angular variable fi in the center of the line {x = 0) for different 
values of the optical thickness tq. Comparing Figs. 2 and 3 we see that the resonance polarization is lower a bit 
than the Thomson one, that is quite expectable. Fig. 4 shows the dependency of P on dimensionless frequency 
X inside the line for different values of tq. We can see that with increasing of tq the line becomes wider, which is 
natural since the number of scatterings rises. The same is true for intensity. The dependency of polarization on 
for To = 0.5 and different frequencies is shown in the Fig. 5a. We can see that this dependency is very sharp in 
the wings of the line, which is explained by small optical thickness there. Finally, the dependence of P on x for 
different values of /i is shown in the Fig. 5b. With the decrease of /i the line becomes wider since the optical path 
along the line of sight increases. 

For To <C 1 we can obtain the following expansions for the components of the vector source function: 

^ ^^=(ro - 2t), 
^ V Zn 



X W ( 1 

X 

+ 2 



To(ro - 2r)(C - 1/6) - In — + (tq - r)^ In ■ 
r(ro-T)ln^^l^ + ^(To-2r) 

T 2 



(28) 



Substituting these expansions into eqs. H26|l and (|27|l we obtain the following expressions for the CMB intensity 
change and polarization profiles (for the observer reference frame): 



r{x,iJ,To) 



(A/2)ro/V2^, ^i < To^{x), 
(/)(x)to, /I > To(/)(a;), 



(29) 



P(x,^,To) 



3 XW_ 



■ 1 

.73 

-(f> (x)' 



64 TT 



C 



In 



To 



In 



^ < To(f>ix), 



To 



^ > To(j>{x). 



(30) 



At low values of tq the last formula is in a good agreement with the results of numerical calculations. For 
example, the polarization dependency on tq is quite sharp (~ Tq) for /i 3> tq(I>{x), but is more slow (~ Tq) for 
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Figure 5: Profiles of polarization in a spectral line: a) for different frequencies x (marked near the curves), b) for 
different values of ^ (marked near the curves). 

l-j, <C TQ(f)(x). For approximate estimation of polarization the accuracy of eq. H3U|I is quite enough. For example, at 
To = 0.01 and /i > 0.1 its error is less than 3.5%. 

4 Conclusion 

In this work we investigate the possible observational display of Thomson scattering on free electrons and of 
resonance scattering in a spectral line of cosmic microwave background radiation in the flat moving layers from 
"dark ages" epoch (10 < z < 1000). The formulae for appearing anisotropy and polarization of CMB are obtained 
and numerical calculations are made. We show that for characteristic parameters according to recent models 
of substance evolution at that period the values of intensity fluctuations and polarization of cosmic microwave 
background may be in the following ranges: A///o = 10"'' 4- 10"^, Q/I = 10^^ 10^^. The effects are observable 
and their discover would be a sufficient forward step in investigations of substance evolution in the pre-galactic 
epoch. 

We do not consider a transfer of the calculated intensity and polarization changes through a homogeneous 
universe from z w 150 to z = 0. The possible distortions caused by such a transfer are defined by an optical 
thickness of the universe due to Thomson scattering. According to the recent measurements in the BOOMERANG 
and MAXIMA experiments (see US]) the Thomson optical thickness of the universe between z — and 

z = 1100 is evaluated to be less than 0.1. As we conserned by an interval between z — and z « 150 the 
distortions of the calculated fluctuations will not exceed probably a few percents which is undoubtedly less than 
an uncertainty of the protoobject models. 
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